1/f noise in variable range hopping conduction. 
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1/f noise induced by traps consisting of donors with no neighbors with close energies in their 
vicinity is studied. Such donors slowly exchange electrons with the rest of conducting media. It is 
shown that in the variable range hopping regime 1/f noise exponentially grows with the decreasing 
temperature. At high temperatures, when the variable range hopping crosses over to the nearest 
neighbor one, we predict a very weak temperature dependence in spite of the activation dependence 
of the conductivity. 
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I. INTRODUCTION 



Low frequency 1// noise was found in many conduct- 
ing materials with the wide variety of transport mecha- 
nisms. Substantial attention was devoted to 1/f noise in 
the hopping conduction §,|,||,|||,||Jl|,|ll| . One of 
the reasons is that 1 // noise limits performance of impor- 
tant devices working in the hopping regime, for example, 
thermistors for X-ray detection and other astrophysical 
applications || . Understanding of the nature of 1/f noise 
can also shed light on the role of electron-electron inter- 
actions in the hopping transport. 

In a lightly doped semiconductor at low temperatures 
electrons are localized on donors and the conductivity is 
due to hopping between donors. (For concreteness we 
talk about n-type moderately compensated semiconduc- 
tor with concentration of donors Njj-) There are two well 
known regimes of hopping transport. At relatively high 
temperatures electrons use almost all the donors for hop- 
ping. This regime is called the nearest neighbor hopping 
(NNH). In this case, the energy scatter of donors due to 
random potentials of charged donors and acceptors plays 
only a secondary role, leading to a fixed activation energy 
of conductivity. At low enough temperatures activation 
required for the use of the majority of donors becomes 
very costly and only donors within a narrow band of en- 
ergies around the Fermi level participate in the conduc- 
tivity. With decreasing temperature the width of this 
band shrinks and hops become longer. This explains the 
name "variable range hopping" (VRH) of the regime. 

Most of recent experiments deal with VRH and focus 
on temperature dependence of 1/f noise ||QJ|||. Re- 
sults seem to be quite controversial. For example, Lee 
|1 found that in silicon 1/f noise amplitude decreases 
when temperature goes down, while McCammon Q ob- 
served exponential increase of 1/f noise with the decreas- 
ing temperature. 

On the other hand, the original theory of 1/f noise 
in the hopping transport QJ^] deals only with NNH and 
only with the case of relatively high temperatures when 
NNH conductivity is temperature independent. This the- 
ory is based on the idea of traps provided by rare isolated 
donors. Such a donor traps an electron from donors of the 
transport paths (conducting media) and releases it back, 



both with characteristic times, which are much larger 
than times of hops determining VRH conductivity. The 
resulting modulation of number of "conducting" carri- 
ers at low frequency lo leads to the spectral density of 
current noise 7^ with a behavior close to 1/u (or 1/f). 
This idea is similar to the McWorter's explanation of 1/f 
noise in MOSFETs by electron exchange between two- 
dimensional electron gas and traps in the oxide, which 
have an extremely wide spectrum of times of tunneling 
to the interface |12J[l|]. The difference is that in the 
MOSFET case traps are separated from the conducting 
two-dimensional gas by the interface, while in the case 
of the hopping transport traps are located inside pores 
scattered in the body of the conducting media of donors. 
This paper generalizes the idea of Refs. [QJ|] to the VRH 
conductivity (Sec. II) and to the temperature dependent 
NNH (Sec. III). We again identify those donors, which 
can work as traps with extremely wide spectrum of very 
large relaxation times. Then we calculate the probability 
of a trap, which now becomes an exponentially decreas- 
ing function of temperature. Integrating contributions 
of all traps to the noise we arrive at the conclusion that 
in the VRH case the spectral density of current noise 7^ 
approximately obeys the Hooge law, Eq. (||), with the 
coefficient a, which is given by Eq. (Q) and exponen- 
tially decreases with temperature. On the other hand, in 
the case of NNH conductivity, which grows exponentially 
with temperature with a constant activation energy, the 
Hooge's a is almost temperature independent. Thus, ex- 
ponential temperature dependence of conductivity does 
not automatically lead to exponential dependence of 1/f 



Sec. IV is devoted to the discussion of the relation- 
ship between of spectral densities of fluctuations of the 
current and the electron concentration, which is used in 
Sec. II and Sec. Ill to calculate the noise spectral density. 
In Sec. V. we compare our theory with two previously 
published theoretical works ||[l0| on 1/f noise in VRH 
and with available experiments |^||,[3||(|. We also make a 
comment about applicability of our theory to MOSFETs. 



1 



II. VARIABLE RANGE HOPPING 

Due to the Coulomb interaction of localized electrons, 
the density of states of donor states, g(e), of a three- 
dimensional lightly doped semiconductor is known [[l4|jl5) 
to have the Coulomb gap near the Fermi level, 



g(e) = (3/^e 



(1) 



Here e is the energy of a localized electron on a donor 
calculated from the Fermi level, k is the dielectric con- 
stant of the semiconductor and e is the proton charge. 
At low temperatures VRH conductivity <r(T) obeys the 
Efros-Shklovskii (ES) law MM 



a{T) = a exp[-(T /T) 



l/2l 



(2) 



where To = Ce 2 /ksna, a is the localization length (Bohr 
radius) of an electron on a donor, and C ~ 2.7. 

The ES conductivity is determined only by donors with 
energies in the band 5 — fcs(T)T) 1 / 2 around the Fermi 
level (5-band). Similarly to Ref. (7]|| we are interested 
in a trap, i. e. in such a rare donor of the same band 
<5, which anomalously slowly exchanges an electron with 
the majority of donors of this conducting band. The rate 
of electron hops between two donors i and j is [fL5| 



vq exp 



2r< 



k R T 



(3) 



where Si is the energy level of the donor i, Sij — 
(\ei\ + \ej\ + \ei — Sj\)/2. We want to find a trap with 
relaxation time larger than v^ 1 , where v is so small 
that ]h{vq/v) » (Xb/T) 1 / 2 . Such a trap is an isolated 
donor, for which all rates of transitions to the neighbor- 
ing donors of the band of energies 5 are smaller than v. 
Let us consider a donor i located in the center of a sphere 
of radius R{y) = (a/2) \a{yo/v), which contains no other 
donors j within <5-band. Then there are no direct hops 
from our donor to the conducting <5-band of donors. But 
this is not enough to have a trap we are looking for, be- 
cause an electron can still escape from donor i to the con- 
ducting 5-band through an intermediate donor j, which 
has energy Ej much larger than 6, but is located closer 
to i, at the distance r.y < R{v). In order to preserve 
the trap all such " dangerous" intermediate donors j near 
donor i should be eliminated or, in other words, all the 
remaining ones should satisfy inequality 



, r > ln — • 

kbJ- v 



(4) 



This inequality means that the energy Ej should be away 
from the energy band of the width A = fc^Tln^/V) 3> 
5 around the Fermi level (A-band). The concentra- 
tion of donors in A-band is of the order of N(A) ~ 
g (A) A = (cA/e 2 ) 3 . Therefore, the average number 
of A-band donors in the sphere of radius R{v) equals 



M = (47r/3)i? 3 (j/)7V(A). Probability, W(v), to find a 
sphere with radius R{u) empty from such donors has a 
form exp(— M) or 



W{v) = cxp ( -B(T/T ) 3 In 6 — 



(5) 



where B is a numerical factor. For a donor with the re- 
laxation rate v equal the typical hopping frequency of 
conducting electrons, v c ~ z/ exp[— (To/T) 1 / 2 ], we get 
W{v) ~ 1, confirming that this is a typical donor. Traps 
have v <C v c and, therefore, exponentially small probabil- 
ities. To find the coefficient B we calculated probability 
of a pore in four dimensional space of three coordinates 
and energy. This calculation gives B ~ 0.3. Note that 
probability W(v) decreases with increasing temperature 
because at high temperatures the volume of the four di- 
mensional pore from which one has to eliminate donors 
becomes larger. 

The spectral density of fluctuations of the concentra- 
tion of "conducting" electrons, n 2 , can be evaluated by 
integration over all traps |L3] 



N(S) 
V 



dW Av 

dv LU 2 + V 



-dv. 



(6) 



Here N(6) is the concentration of donors in the 5-band 
around the Fermi level and V is the sample volume. 

According to the standard approximation is that rela- 
tive fluctuations of the concentration of conducting elec- 
trons lead to fluctuations of conductivity, i.e. 



P N 2 (S) 



(7) 



Here I 2 is the spectral density of current fluctuations and 
/ is the average current through the sample. We will dis- 
cuss justification of Eq. (0) in Sec. IV. Using Eq. (||) 
we can now estimate I 2 /! 2 - Let us start from the most 
interesting case when the exponential factor dW/ dv in 
the integrand of Eq. (^|) is a weaker function of z/, than 
v '/{ui 2 + v 2 ) and the integral is determined by v ~ u>. 
Then one can take exponential term of Eq. (](]) out of the 
integral at the frequency v = u>. This leads to the Hoogc 
law lO 



I 2 ujN d V 1 



(8) 



where NjjV is the total number of donors and a is the 
Hooge's coefficient. We obtain 



a(w,T) oc exp[-B(T/T ) 3 In 6 — ]. 



(9) 



Here as everywhere below we are concentrated on ex- 
ponential temperature dependence of noise and are not 
trying to calculate prefactor. Using Eq. ( p"2| ) one can 
express a(u,T) through <r(T): 
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-^--Kra)'- (10) 

This simple relationship can be convenient for a direct 
experimental verification. 

The above derivation of ce(ui,T) is justified when the 
dependence of « on w is weaker than l/u>, i. e. when 
u > ^ exp[-(T /T) 3 / 5 ]. On the other hand, our ap- 
proach of traps is good only for frequencies u> smaller 
than the frequency of transport hops v c , i. e. for 
\n{v Q /uj) > (To/T) 1 / 2 . Thus, the Hooge law is valid 
in the range of frequencies 



exp[-(T ( 



o/T) 1 / 2 ] » u/vq » exp[-(T /T) 3 / 5 ] 



(11) 



In spite of close powers in the exponents of lower and 
upper limit this range can be quite broad. For example, 
even at a moderate (To/T) 1 / 2 = 10 this range covers Z.b 
decades, while at the very large (To/T) 1 / 2 = 20 it reaches 
7 decades. 

One can approximately write I 2 /I 2 oc 1/w 7 , where 
7=1- 6B(T/T ) 3 lnVoM. Thus, according to this 
theory 7 < 1. In the range (|ll|), the index 7 is close to 
unity and approaches 1/2 at the low frequency limit of 
applicability of Eq. (||). 

For a fixed uj in the frequency range ((H)), Eq. (||) 
contains the main result of this paper - an exponential 
growth of the noise power with the decreasing tempera- 
ture. The physics of this dependence is clear: the lower 
the temperature, the smaller is the width of the band 
of energies A, where one has to eliminate 'dangerous" 
donors in order to get a trap for a given frequency, the 
smaller is the number of donors to be eliminated, the 
larger is the probability to find such a trap. 

According to Eqs. (]?]) and (0) at smaller frequencies 
oj <C Vmin = va exp[— (To/T) 3 /^], the spectral density of 
noise becomes substantially weaker than 1/ / and eventu- 
ally saturates at the level ~ exp[(To/T) 3 / 5 ]. To calculate 
behavior of 1// noise in all range u> < v c more accurately 
one should similarly to add the contribution of traps 
made of large finite clusters of donors and calculate the 
integral of Eq. (||) numerically. 

We derived Eq. (g) for a semiconductor in which the 
concentration Nd is substantially smaller than the criti- 
cal concentration of the metal- insulator transition, Nmi- 
In this case, we can use Eq. (|l|) for the density of states 
g(s). In experiments VRH is typically measured at 
0.25Nmi < Nd < Nmi because of a huge resistance 
of lighter doped samples. We believe that Eqs. (||) and 
@ still work in the range 0.25N MI < N D < 0.5N MI . 

We do not exactly know how to gener alizc Eq. (g) 
for the critical vicinity of the metal-insulator transition, 
where Nmi ~ Nd <C Nd- An additional source of the 
traps may be related to the spacial fluctuations of diverg- 
ing at the transition localization length £. Anomalously 
small £ in a vicinity of some donor can create a trap. A 
simple conjecture is that Eqs. (||) and (jl^) remain valid 



at Nmi—Nd < N D with the same T = 2.7e 2 /fcs«;f as in 
the ES law near the transition. Anyway an experimental 
study of this range of concentration can shed some light 
on the metal-insulator transition. 

Until now we have been dealing with VRH conductivity 
in a doped crystalline moderately compensated semicon- 
ductor, where conductivity obeys Eq. ([l2]). In systems 
with a weaker role of the electron-electron interaction 
such as amorphous semiconductors or two-dimensional 
systems screened by a parallel metallic gate, VRH can 
take place in the larger band of energies than the width 
of the Coulomb gap and, therefore, one can observe Mott 
law 



a(T) =a Q exp[-(T d /T) 



X/(d+l)l 



(12) 



where d = 2,3 is the dimensionality of the system, T<j — 
(3d/{ga d ), g is the bare (unperturbed by the Coulomb in- 
teraction) d-dimensiona! density of states in the vicinity 
of the Fermi level and (3d are numerical constants [ JL5| . 
Repeating the derivation of Eq. (g) we get the Hooge law 
with 



a(uj,T) ~ exp[-B d (T/T d )\n a+1 (v /u)] , 



lna(o;, T) = — Bd 



( ln(i>o/(j) 
\ln(a Q /a) 



d+l 



(13) 



(14) 



where Bd are numerical coefficients. The physics of the 
origin of the exponential temperature dependence is ex- 
actly the same as above: the higher the temperature, the 
more difficult is to find a trap for a given small frequency 



III. TEMPERATURE DEPENDENT NEAREST 
NEIGHBOR HOPPING 

Eqs. (g) is valid when temperature is so small that 
the energy width A of the band, where all "dangerous" 
donors are eliminated, is smaller than the width of the 
donor band, A = e 2 Np 3 /n, i. e. at fceTln^o/^) <C 
e 2 Np 3 /k. In the opposite case, A 3> A, or 



T>T = 



e 2 N D /3 



fcsKln(^o/^) ' 
the concentration of donors in the A-band 

N(A)= / g(e)de = N D [l-(A/A) 3 } 

J-A 



(15) 



(16) 



is close to Nd and only weakly depends on A. Here 
we used the fact that as shown in Chapter 3 of Ref. 
pi the large energy tails of the density of states of 
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the classical donor band behave as g(e) ~ NdA 3 /\e\ 4 
at |e| >• A. Now we can use Eq. (||) with probabil- 
ity W{v) = exp[-( 7 r/6)A^(A)a 3 ln 3 (t'o/^)] which can be 
rewritten as 

W(u) = exp[-(7r/6)N D a 3 ln 3 ( V)] exp(T c /T) 3 . (17) 

Here T c ~ e 2 aN^ /3 /k B K is the critical temperature of 
the transition from NNH to VRH In the range of 

frequencies 

expHTo/T) 1 / 2 ] > u/vo » exp[-{N D a 3 )-^ 2 ] , (18) 

where dW/dv is weaker function of v than 4i//(cj 2 + v 2 ). 
we again obtain the Hooge law with 

a(u;,T) =exph(V6)iVDa 3 ln 3 (^/ W )]exp(T c /T) 3 . 

(19) 

Eq. @ is valid both for NNH (at T > T c ) and for VRH 
(at Ti < T < T c ). First exponential factor of Eq. |[|) 
originates from the probability of a pore free of all donors. 
This factor is temperature independent and coincides 
with the value of a obtained for the high temperature 
limit of NNH (7||] . Second, temperature dependent expo- 
nential factor of the right hand side of Eq. ( |l9| ) provides 
only a relatively small correction to lna(w,T), which 
grows as 1/T 3 with the decreasing temperature. In NNH 
range, when T 3> T c , the relative temperature dependent 
correction to a{u>,T) is much smaller than unity. At the 
transition point from NNH to VRH, i. e. at T ~ T c , this 
correction reaches 100%. At T <C T c it becomes exponen- 
tially large and provides the smooth crossover to Eq. (||) , 
when T approaches T\ — e 2 N D { 3 /fes/tln^o/i;). How- 
ever, as we already mentioned above, because of huge 
resistances an experimental observation of VRH is not 
possible, when N^a 3 <C 1. Therefore, inequality ( |l8| ) is 
very difficult to realize and the use of Eq. (|l^) for VRH 
regime is limited. But Eq. @ is definitely important 
for the temperature dependence of noise of NNH, which 
takes place at N^a 3 -C 1. 

To conclude this section we would like to repeat that 
for NNH temperature dependence of a(u>, T) is negligible 
in all range of temperatures A 3> T 3> T c , where con- 
ductivity itself grows with decreasing temperature with 
activation energy of the order of A [|15| . 

IV. RELATIONSHIP BETWEEN 
FLUCTUATIONS OF THE CURRENT AND OF 
THE CONCENTRATION OF CONDUCTING 
ELECTRONS 

Now we would like to return to Eq. (^), which estab- 
lishes relation between fluctuations of concentration of 
conducting electrons (which are not trapped) n and fluc- 
tuations of current and discuss its validity. In MOSFETs, 



where most of electrons are free, fluctuations of the con- 
centration of free electrons obviously lead to proportional 
fluctuations of the conductivity [ ^2|JT^ ] and at a fixed elec- 
tric field 11/ 1 2 = nl/n 2 . 

For NNH such proportionality between fluctuations of 
/ and n is obvious for strongly compensated samples, 
where electrons play the role of carriers, or for a weakly 
compensated semiconductor, where the role of carriers is 
played by a small concentration of holes (empty donors). 
However, even in this case there could be such an in- 
termediate compensation ratio K = Na/Nd, where the 
coefficient in linear proportionality between fluctuations 
of / and n vanishes. (Here Na is the concentration of ac- 
ceptors). In VRH conductivity connection between fluc- 
tuations of I and n is even less trivial. 

Let us consider relation of I and n for a simpler ex- 
ample when concentration, n, of electrons in a two- 
dimensional hopping system is varied due to a change 
of the gate voltage. Any variation of n leads to a shift 
of the Fermi level. The Coulomb gap, however, moves 
together with the Fermi level and, in the first approxi- 
mation, it remains unchanged at small energies 8, which 
are important for VRH. It is not obvious then whether 
Eq. @j is valid. 

In the first subsection of this section we show that in 
the generic case of asymmetric with respect of Fermi level 
bare (disorder related) density of states one can justify 
Eq. (0), because in this case there is a nonzero derivative 
da/dn for a sample of infinite size. In the second sub- 
section we introduce more general mechanism of fluctu- 
ation of conductivity, which exists even in the case when 
da/dn = for infinite sample, but still can be evaluated 
in the way it was done in Sees. II and III. 



A. Asymmetry of the density of states 

In a doped n-type semiconductor at a generic compen- 
sation ratio K the density of states at large e is asymmet- 
ric with respect to the Fermi level (see plots in Chapter 
14 of Ref. fL5fl). In spite of tendency of Coulomb gap 
to be symmetric, the density of states, 17(e), at e ~ 6, 
is sensitive to behavior of the density of states at larger 
energies. This happens because the Coulomb gap should 
crossover in its " shoulders" to the bare density of states of 
the classical impurity band |H| . Using the self-consistent 
equation for the one-electron density of states derived by 
Efros [fl6| one can show that the absolute value of nonuni- 
versal (depending on position of Fermi level) correction to 
the density states at small s is proportional to K 3 e 3 /e 6 A, 
where A ~ e 2 A^ /3 / k is the energy width of the classi- 
cal donor band. In other words, the relative asymmetric 
correction to Eq. ([j]) at e ~ 5 is proportional to 8/ A. 

Our final goal is to evaluate fluctuations of the hop- 
ping conductivity. In VRH regime the conductivity 
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fluctuates because g(S) fluctuates following the Fermi 
level. It is easy to show that relative fluctuations of 
the conductivity are larger than fluctuations of g(5) by 
(Tq/T) 1 /' 2 ~ TqJ 8. Thus, as in a free electron gas, the 
fluctuations of conductivity are proportional to fluctua- 
tions of the concentration of conducting electrons with 
the temperature independent coefficient of proportional- 
ity (6/A)(T /S) = T / A. 

Now we can return to slow traps immersed in the con- 
ducting media, which work similarly to a gate. While the 
number of trapped electrons is slowly changing with time, 
the conducting media arrives at the quasi-equilibrium 
at much shorter time scale and develops the quasi- 
equilibrium Fermi level tracking fluctuations of the elec- 
tron concentration. In turn, the position of the quasi- 
Fermi level affects the density of states relevant to the 
conductivity and, therefore, the hopping conductivity it- 
self. 

Without Coulomb interactions fluctuations of the 
concentration are uniform over the conducting media. 
Coulomb interactions makes these fluctuations localized, 
because of screening of trapped electrons. When an elec- 
tron is trapped, a hole remaining in the conducting media 
is localized close to the trap. In other words, fluctuations 
of the concentration of conducting electrons do not prop- 
agate to the bulk of the conducting media. In this sense, 
traps are equivalent to small gates separated from each 
other by large distances. They locally perturb the con- 
centration of conducting electrons. This perturbation, 
of course, changes conductivity only near traps. Now we 
have to discuss these changes and the way how they affect 
the total current through the sample. 

We start from an estimate of the screening radius of 
a spherical trap with radius R = (a/2)hi(vo/v). In 
other words, we want to find the radius of the concentric 
sphere, where the hole left behind by trapped electron 
is localized. In the Coulomb gap, the density of states 
at the Fermi level is very small at small temperatures. 
One can find substituting ksT for e into Eq. (H) so that 
linear screening Debye radius R s which describes screen- 
ing of very small charge is very large, R s ~ cl(Tq/T). 
On the other hand, one electron charge — e of the trap 
is large enough to produce additional nonlinear screen- 
ing. Nonlinear screening radius Rn can be found from 
the condition that potential of trapped electron —e/kR 
can attract a positive charge to one of donors inside the 
sphere with radius Rn- This condition can be written as 

4-7T re/itR 

Y (r 3 N -R 3 )j 2g(e)de~l, (20) 

This gives Rn — R ~ R- To understand changes of the 
macroscopic conductivity we should compare Rn ~ R 
with the characteristic distance L at which VRH con- 
ductivity self averages, or in other words, becomes uni- 
form. This is the characteristic period of the critical 



subnetwork (percolation infinite cluster) of the Miller- 
Abrahams resistor network, which carries all the cur- 
rent p5| . It is known |H| to be of the order of L = 
(a/2)(T /T)( 1+ ^)/2 ~ a(T /T) - 95 . (Here we used the 
fact that exponent of the correlation length of the per- 
colation theory, = 0.9 ~ 1 |l5|.) The length L is of 
course much larger than the characteristic length of a 
hop rh — (a/2)(T /T) 1 / 2 . It is interesting that L is only 
somewhat smaller than the linear screening radius R s . 

Below we consider two different cases, R > L,R S and 
R < L,R S . In the first case, Rn > R s and, there- 
fore, screening is linear ( Rn plays no role). The hole 
is not bound to a particular donor but freely travels in 
the spherical layer (atmosphere) between R and R + R s . 
In this case, one can use macroscopic approach for the 
role of traps, namely one can say that the concentration 
of electrons (and the local VRH conductivity) slightly 
in varies in spherical atmospheres around all traps. It 
is known that in a slightly inhomogeneous media in the 
first approximation the sample conductivity is equal to 
the average value of the local conductivity. Therefore, the 
change of the sample conductivity reflects fluctuations of 
the total number of conducting electrons and Eq. (Q) is 
valid. 

However, the case R > L,R S can be realized only at 
very low frequencies v and we have to consider the sec- 
ond case, R < L, R s . In this case, the simple macro- 
scopic approach based on fluctuations of concentration 
does not work, because traps perturb the critical network 
locally, at the distances much smaller than the scale of 
self-averaging of the VRH conductivity. We show below 
that for R < L,R S l/f noise exists even for the per- 
fect symmetry of the density of states with respect of the 
Fermi level. 

B. Effect of the trap potentials 

In the case R < L, R s nonlinear screening is impor- 
tant because Rn < R s - Therefore, the hole is bound 
to the trapped electron and is localized on a particular 
donor. In other words, electron travels forth and back be- 
tween this donor and the trap creating the dipole poten- 
tial 4>(R,t), which modulates conductances G,j of neigh- 
boring resistors of the Miller- Abrahams resistor network 
p5| . In order to calculate the effect of this potential 
let us devide the sample in cubic blocks with linear size 
L each and consider for simplicity such frequencies u, 
when R < L, but is still comparable to L, say R ~ L/5. 
Only very small fraction of these cubes of the order of 
(dW/dv)v has one trap inside. When an electron is cap- 
tured or released, the dipole potential <fr(R, t) changes 
energy of critical resistors which derermine conductance 
of two halfs of the block on both sides of the trap by 
energy e<j>(R,t) = e 2 / kR ~ ksT. Positive and nega- 
tive potentials of the dipole applied to two different halfs 
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of the cube, produce different effects because two halfs 
have different random conductivities (conductivity self 
averages only at a distance larger than L). As a result 
the conductance of such cube changes by 100%. Thus, 
due to the trapping and detrapping of an electron the 
conductance of the cube with the trap inside fluctuates 
roughly speaking by factor two. Macrosopic conductivity 
then fluctuates according to Eqs. (||) and 



V. COMPARISON WITH OTHER 
THEORETICAL MODELS AND EXPERIMENT 

Let us start from two theoretical papers, which aim at 
the theory of 1/ f noise in the VRH conduction. Kozub ||] 
suggested seemingly different mechanism related to pairs 
of donors slowly exchanging electron with each other and 
modulating transport on current paths by their poten- 
tial. To produce effect at a very small frequency, two 
donors should be far from each other. But they also 
should not be "shortened" by a chain of other donors 
providing faster exchange between them. In other words, 
they should be isolated from the conducting media. One 
could think that each of two donor should have a pore 
similar to ones discussed in our paper around it. In this 
case, for a given frequency a pore for a "modulator" pair 
would have much smaller probability than spherical traps 
we study and "modulators" would produce a negligible 
noise on the background of the noise created by traps. 

Actually, a pair is not "shortened" even if only one 
of its donors has a spherical pore around it or, in other 
words, this donor is exactly a trap we study above. The 
second donor is the donor of the (5-band closest to the 
pore. So one can say that we presented here a quantita- 
tive development of the Kozub's qualitative idea. Indeed, 
the language of the second half of Sec. IV is very close 
to that of Ref. §. 

In another paper, Shtengel and Yu JIo| apparently 
deal with the same mechanism as discussed in our pa- 
per. They did not study analytical asymptotic dependen- 
cies, but evaluated expression for I 2 /I 2 numerically and 
obtained results dramatically different from ours. They 
found stronger than I// growth of I 2 /I 2 at small fre- 
quencies and a weak growth of the noise with tempera- 
ture. I have no explanation for their numerical results. 

Let us now compare our predictions with experiments 
on crystalline semiconductors concentrating mainly on 
the temperature dependence. In the range of VRH expo- 
nential growth of 1// noise with decreasing temperature 
following from Eq. (|^) agrees qualitatively with 1// noise 
measurements for implantation doped silicon samples || , 
but disagrees with results on bulk silicon crystals H . On 
the other hand, in agreement with Eq. (19), 1// noise is 
almost temperature independent for NNH in germanium 
H . I // noise in the wide range of temperatures including 
the transition between VRH and NNH was also studied in 



two-dimensional n-GaAs channels || . In agreement with 
our theory it was observed that I 2 /I 2 sharply decreases 
with increasing temperature in VRH range, while in NNH 
range the temperature dependence becomes weak. 

Frequency dependencies are more problematic for the 
theory, because deviations in the direction of smaller 
power or saturation at small frequencies were not seen in 
experimental works. This can be explained for samples 
with very large resistances. When (Tb/T) 1 / 2 is very large 
the frequency of saturation can become too small to be 
measured. Indeed, at v = I0 12 Hz and (Tq/T) 1 / 2 = 20, 
we find that v c = 2.5 10 3 Hz and v m in = 2.5 I0 _4 Hz. On 
the other hand, the absence of a tendency to saturation 
in moderately resistive samples with (To/T) 1 / 2 ~ 10 may 
signal that the actual mechanism of I// noise is different 
from the one considered here. Indeed, there are many dif- 
ferent pseudoground states separated by large potential 
barriers from each other in the Coulomb glass formed by 
a frozen distribution of interacting electrons on random 
donors Jl^,|Tl|. Different pseudoground states (valleys) 
can have different conductivities and random walk of the 
system through a sequence of these states may lead to 
1// noise |n]] . Logically it is possible that this noise can 
be larger than the one we estimated, but nobody was 
able to evaluate it. 

Until now we have dealt with the hopping transport, 
where both roles of conducting media and traps are 
played by localized states of impurities. There are also 
systems where current is carried by free electrons, while 
localized states play the role of traps only. The best 
example of such a system is a MOSFET, where a free 
two-dimensional gas exchanges electrons with localized 
states (traps) of the oxide. McWorter's theory |l2j as- 
sumes that density of states of traps is so small that any 
localized state of the oxide, which is close to the interface 
and close in energy to the Fermi level plays the role of 
the trap. These traps are assumed to be uniformly dis- 
tributed with respect to the distance from the interface 
and, thus, lead to I// noise. 

McWorter's theory [[l2| neglects the possibility that 
near a trap, which exchanges electrons with the conduct- 
ing channel, other localized states can be located. They 
can provide alternative hopping paths with exponentially 
smaller exchange time and "shunt" the trap. However, 
for small enough frequencies and a large enough den- 
sity of localized states in oxide these paths always exist. 
Then most effective traps consist of a localized state sur- 
rounded by a spherical pore, which touches the interface 
of the channel. As in the case of VRH, all localized states 
within energy band A should be excluded from the pore, 
while localized states with higher energies can stay there. 
Therefore, one can say that the pore is arranged in the 
four-dimensional coordinate-energy space. 

The theory of this paper can be directly used to cal- 
culate the probability of such traps and therefore the 
frequency dependence of spectral density of noise. We 
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assume that in the gap of the oxide the bare density of 
states, g is constant. The Coulomb interaction forms 
the Coulomb gap of the width Eqg ~ e 3 g 1 ^ 2 /n 3 ^ 2 in 
the vicinity of the Fermi level. One has to compare the 
energy "width" of the trap A and E C g- If A < E C g 
one can use Eq. (^) On the other hand, when A > Ecg 
Eq. (0) crosses over to Eq. ([l3|). 

In conclusion, this paper explores the role of isolated 
donors as traps which modulate the concentration of con- 
ducting electrons and create noise of the hopping conduc- 
tivity. It is shown that in the wide range of frequencies 
traps lead to approximately 1/ / behavior of noise, with 
the crossover to saturation at extremely small frequen- 
cies. In the nearest neighbor hopping range, in spite of 
the activation temperature dependency of the conductiv- 
ity, the relative intensity of 1// noise is almost tempera- 
ture independent (see Eqs. (||) and (p^)). In the variable 
range hopping regime the intensity of 1// noise grows 
exponentially with the decreasing temperature according 
to Eqs. (|) and (§). 
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M. Weissman for the critical comments on original ver- 
sion of this paper, which led to understanding of the role 
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